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Cosmic-ray scattering on magnetic turbulence leads to spatial diffusive propagation; if the scatter¬ 
ing medium is moving, this will inevitably also cause changes in the momentum of the particles, 
so-called diffusive reacceleration. This can be described as diffusion in momentum space. Diffu¬ 
sive reacceleration has often been invoked to explain the peak observed in secondary-to-primary 
ratios at a few GeV, in particular Boron-to-Carbon. This avoids the necessity to postulate an 
ad-hoc break in the spatial diffusive coefficient, and has become almost a standard in modelling 
cosmic-ray spectra. However, at the levels invoked, the process implies a significant input of 
energy from the interstellar medium into cosmic rays, so that in such models interstellar space 
competes with the usual accelerators like supernova remnants. The questions arise: is reacceler¬ 
ation really occurring at the high level required to explain secondary-to-primary ratios? and are 
the energy requirements physically plausible? We address this issue using both analytical and 
numerical models of cosmic-ray propagation. 
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1. Introduction 

In most discussions of cosmic ray propagation it is assumed that the cosmic rays are produced 
in discrete sources, generally taken to be supernova remnants, although there is increasing interest 
in and evidence for possible additional classes of sources, see e.g. the presentations in the recent 
conference [p, and then diffuse through the interstellar medium and an extended cosmic ray halo 
before escaping from the Galaxy. This diffusion is not just in space, but also in momentum, if, 
as is inevitably the case, the magnetic fields that scatter the particles are not purely static but are 
in random motion with characteristic velocities of order the Alfven speed. This leads to a certain 
amount of second order Fermi acceleration during propagation, an effect which for historical rea¬ 
sons is usually referred to as “reacceleration”. There has been a long and inconclusive debate as 
to the importance of reacceleration, but there is no doubt that the idea is attractive for a variety of 
reasons, and physically the effect must occur at some level. 

One of the strongest arguments in favour of reacceleration is that it allows a rather natural fit to 
the low energy Boron to Carbon data. This was persuasively argued by |§] who showed that with 
reacceleration the B/C data could be well fit with a single power-law dependence of the escape 
path-length as a function of rigidity, whereas in a leaky-box model one had to suppose an unnatural 
decrease in the path-length at low energies. The fact that the inferred power-law dependence of the 
path-length also agreed with that predicted for Kolmogorov turbulence was seen as another positive 
feature of the model (although whether interstellar magnetic turbulence should be described by a 
Kolmogorov spectrum is another question). 

This reacceleration by interstellar turbulence is a potentially significant contribution to the to¬ 
tal Galactic cosmic ray source power (and of course also a potentially significant damping term 
for interstellar turbulence). This was qualitatively discussed in [^] on the basis of analytic esti¬ 
mates. In this paper we seek to make quantitative estimates of the reacceleration power using the 
well-established cosmic ray propagation code Galprop as first described in [Q]. Galprop has 
become the de facto reference model for cosmic ray propagation studies in our Galaxy and is thus 
an appropriate choice, but any other propagation model incorporating reacceleration should give 
very similar results. The power transferred to the cosmic rays comes of course at the expense of the 
general turbulence in the interstellar medium so we also evaluate a damping process for interstellar 
turbulence. Finally we note that the whole question of propagation models has acquired new rele¬ 
vance in the context of recent searches for possible signals of dark matter in cosmic rays, see e.g. 
[^] and references therein. 

2. Reacceleration power 

This section briefly reviews and expands on fhe main resulfs of [||]. The basic equation under¬ 
pinning all diffusion models of cosmic ray fransporf, going back af leasf fo fhe classic monograph 


of [|] i 


IS 



( 2 . 1 ) 


Here f{x,p,t) is fhe isofropic part of fhe phase space densify of a given species as a function 
of position X, scalar momenfum p and time t and Q{x,p,t) is a source term representing the initial 
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production and acceleration of cosmic rays. The spatial and momentum space diffusion coefficients 
are D^x and Dpp respectively and will in general be functions of position, momentum, and time. 
The ellipsis denotes the many other terms that should be included, for example those describing 
ionisation energy loss and nuclear interactions, but which are irrelevant for this discussion. 

On heuristic grounds ||3|] the two diffusion coefficients are related by 


Dxx^pp 


-yvi 


Detailed calculations |^, see also equations 9.38 and 9.39 in [^, show that 

l-At^ \ /(l-At2)v+v_' 


DxxDpp = py 


Vi +V_ 


v+ + v_ 


( 2 . 2 ) 


(2.3) 


where the angle brackets denote an average over an isotropic pitch angle distribution p being the 
pitch angle cosine. 



(2.4) 


and v± is the rate of pitch-angle scattering off forward and backward propagating Alfven waves. If 
the scattering rates are taken to be constant and equal we recover the naive heuristic result, but in 
general the scattering rates should be functions of the pitch angle and related to the wave spectrum 
through the gyro-resonant condition. This introduces a weak dependence on the shape of the wave 
spectrum, but not on its amplitude. Specifically, if we assume a power-law wave spectrum with 
W (k) oc where jW{k)dk is the total energy content of the waves and k is the wave-number 
(a = 5/3 being a Kolmogorov spectrum for example). 


^reslf^ib (^res) 


where Vg = p/eBo is the particle gyroradius, Bq the background magnetic field, and the resonant 
condition is 

prgk^Qs ^ 1 . ( 2 - 6 ) 

Considered as a function of pitch angle the scattering rates thus scale as 


v± oc /i" ’ 


(2.7) 


so that (assuming v+ = v_) 

/ 1-At^ \ / (l-At^)v+v_ \ 

\ V+ + v_ / \ v+ + v_ / 

= J f\l-pV-^dp 
4 Jo Jo 

= -0<a<2 

4\2 — a 4 — aJ\a 2-\-aJ 

1 

fl;(4 — a) (4 — a^) 
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Thus more generally 


D„Dpp = p^vl- 


(2.9) 


a(4 — a) (4 — a^) 

which agrees with the simple heuristic estimate for a = 1. 

However if we hold jx fixed and look at the scattering rate as a function of particle momentum 
or gyro-radius, then 

V± oc vr„'^W(kres) vr^g (2.10) 


and thus 


Dxx°^V^( (l-ft^) 


1 


V+ + V- 


p\ {2-a) 

oc V 1 — 

e, 


( 2 . 11 ) 


which is the commonly used form, D oc vR^ or velocity times a power-law in particle rigidity 
(momentum per charge, which is directly proportional to the gyro-radius for fixed magnetic field). 
The power-law index of fhe rigidify dependence of fhe spafial diffusion coefficienl, 5, is fhus related 
fo fhaf of fhe wave specfrum, a, by 


5 = 2 — a, a = 2 — 5 


( 2 . 12 ) 


and, as is easily verified, 

1 1 

fl(4-a)(4-fl2) - 5(4-5)(4-52)- 

In facl fhe denominafor can be written 

a(4 —a)(4 —a^) = 9 — 10(a — 1)^ -|- (a — 1)“^ 


(2.13) 


(2.14) 


which shows fhaf if is symmefric abouf and has a local maximum af a = 1. Thus wifhin fhis simple 
quasilinear fheory of particle scattering we have 


DxxDpp = p^Vl 


5(4-5)(4-52) - 9^^^^ 


(2.15) 


which is essenfially fhe form used by Galprop, see equation (1) in derived from []10|], if fheir 


factor 4/3w is sef fo unify (fhe inclusion of w appears formally incorrecf, buf as if is freafed as a 
consfanf of order unify if makes no difference fo fhe resulfs). If is interesting fo note fhaf for equal 
infensifies of forward and backward propagafing waves fhe naive scattering esfimafe corresponds 
fo fhe minimum possible level of reacceleration, buf of course lower values can be obfained wifh 
anisofropic wave fields. 

For convenience we will from now on wrife 


Dpp = i»p^vl 


1 

Dxy 


(2.16) 


For an isofropic Kolmogorov specfrum wifh 5 = 1/3 we have r? = 81/385 « 0.21 and fhe largesf 
if can be is for 5 = 1 where we have r? = l/9RiO. 11; fhe resulfs of 0] can be obfained by setting 


t? = 


4 

35(4-52)(4-5)w' 


(2.17) 
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The local reacceleration power density is shown in [Q], assuming very reasonable regularity 
conditions on the particle distribution function, to be given by the integral over the particle spec¬ 
trum, 

1 d 

Pr = Anp^f—— [p^Dppv) dp (2.18) 

which can be written, if Dpp is expressed in terms of D^x using the above relation, as 


Pr 



' d\n{v/Dxx) 
d\np 


dp. 


(2.19) 


If we parametrise the spatial diffusion as 



( 2 . 20 ) 


where m is the particle mass and c the speed of light this can be written in the useful and rather 
transparent form 

Pr = Td(4 — 5)—mc^ [47lp^f(—\ dp (2.21) 

Dq J \ me / 

which gives the local reacceleration power density in terms of a simple integral over the spectrum. 
Integrating this over the entire Galaxy then gives the total power transferred to the cosmic ray 
population from Alfvenic turbulence. 


2.1 Analytical estimate 

The obvious first step is to consider a power-law distribution function of the form 


f{p) = fo 



P > Pmin 


( 2 . 22 ) 


above some minimum momentum pmin where observationally y k, 4.% (the distribution function in 
momentum is steeper by 2 than the differential energy spectrum, so this corresponds to N{E) oc 
£■- 2.8 fQj. relativistic particles). However as is easily seen this gives an integral which is dominated 
by the low-energy cut-off. 


Pr 


4np^fo — 


Pmin 


me 


r/2 

t?(4 — 5)—mc^ 

Do 

rJ(4-5)^mP(mc)Vo^q:^ 


1 - 5-7 


dp 

V me J 


(2.23) 


(note that 5 -I- 7 > 4) and diverges as pmin —> 0 

The total reacceleration power is thus dominated by the poorly-know low-energy part of the 
cosmic ray spectrum, as noted in [|^, but we can at least estimate the contribution from particles 
above about 1 GeV/n where we have reliable data. 

The differential energy flux of cosmic rays above 1 GeV is given by the Particle Data Group 
booklet as approximately 

/ E \ 

1(E) = 1.8 X 10"^ - m-^s-i sr-i (GeV)-i (2.24) 

V 1 GeV / 
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corresponding to (for simplicity, and as this is a rough estimate, we consider only protons and take 
the proton rest mass to be 1 GeV) a proton distribution function of 

= - = 6 X 10^^ f—) m-hi-UGeYy^ (2.25) 

c \mc/ 

This implies an energy density in relativistic cosmic ray protons of 

lAKp^mpcdp^iKlO-Hnmc^f^ (£) 

Ri l.OeVcm^^ (2.26) 


and a reacceleration power density for relativistic cosmic rays of 


Pr = T^{4-5)^mc^ [4np^f(—) dp 

Do J \mcJ 

^ ^{4-5) 

= t? 




1 


Do 0.7+ 5 
4-5 / 14 


0.7 + 5 V 30kms 


0.75 eVcm 

2 / 


-3 


Do 


V lO^^cm^s ^ 


-1 


7 X 10^‘®eVcm^^s^' 


(2.27) 


with more from the low-energy part of the spectrum. For a typical CR residence time in the Galaxy 
of at least 10^ yr, reacceleration thus adds at least 20% to the energy content. To proceed further 
we need to turn to numerical methods. 


2.2 Estimate using GALPROP 

We use the numerical CR-propagation package GALPROP\ using parameters which repro¬ 
duce B/C with reacceleration, and also with the same parameters but without the reacceleration 
term. The total CR proton energy content of the Galaxy is computed by integration over momen¬ 
tum and volume. The luminosity is computed assuming a nominal residence time in the Galaxy. 
We base the calculation on the model z04LMS, which has reacceleration and a CR halo height of 
4 kpc l|TT]]. The effective CR residence time is here set to 3 x 10^ yr to reproduce the total proton 
luminosity of 8 x 10'*° erg s^* [jlT]], but the relative values are independent of this assumption. 

For Va = 30kms^* the proton luminosity is 7.8 x lO'^^ergs^* and for Va= 40 km s^* , 8.7 x 
10^*° ergs'*. Without reacceleration, the proton luminosity is 6 x 10“*° erg s'*. These are typical 
parameter values which reproduce B/C in reacceleration models (see also [|^]), so 30-45% of 
energy input is from reacceleration. This is consistent with the analytical estimate above. 


3. Conclusions 

Our estimates of the reacceleration power are necessarily approximate, but suffice to demon¬ 
strate that, for parameters chosen to reproduce B/C, as is often done, the energy input via reac¬ 
celeration is comparable to the total energy, which is usually considered to arise in SNR. SNR 
are then not “the sources of CR” (in the sense of the sole source of their energy) in such models 

'Current version available at https://sourceforge.net/projects/galprop 
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(although of course they can still contribute indirectly as being among the many agents driving 


the bulk turbulence of the interstellar medium []13|])! The questions of whether this is physically 
plausible, and what experimental tests could be made to confirm reacceleration or otherwise, are 
left to a future paper. We note however that the main dissipation channel for interstellar turbulence 
is almost certainly the non-linear cascading to small scales (indeed this idea is implicit in the use 
of Kolmogorov or Kraichnan spectra) and thus the total power requirement for maintaining the in¬ 
terstellar turbulence could be many times the energy going into diffusive reacceleration of cosmic 
rays. It is possible that the cascading simply transfers the energy to the relevant scales for reaccel¬ 


eration and that indeed cosmic ray reacceleration is the dominant dissipation process, see [14], but 
then one can certainly not use power-law wave spectra to fit the B/C data. 

Meanwhile it is clearly desirable to consider alternative explanations for the peak in the energy- 
dependence of B/C, rather than blindly fitting the data with reacceleration; in particular convection 


gives a fairly natural mechanism, and there is evidence for a Galactic wind [15]. 
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